Under consideration for publication in Phys. of Fluids 

Effect of the Earth's Coriolis force on the large-scale circulation of turbulent 

Rayleigh-Benard convection 
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We present measurements of the large-scale circulation (LSC) of turbulent Rayleigh-Benard con- 
vection in water-filled cylindrical samples of heights equal to their diameters. The orientation of 
the LSC had an irregular time dependence, but revealed a net azimuthal rotation with an average 
period of about 3 days for Rayleigh numbers R > 10 10 . On average there was also a tendency 
for the LSC to be aligned with upflow to the west and downflow to the east, even after physically 
rotating the apparatus in the laboratory through various angles. Both of these phenomena could be 
explained as a result of the coupling of the Earth's Coriolis force to the LSC. The rate of azimuthal 
rotation could be calculated from a model of diffusive LSC orientation meandering with a potential 
barrier due to the Coriolis force. The model and the data revealed an additional contribution to the 
potential barrier that could be attributed to the cooling system of the sample top which dominated 
the preferred orientation of the LSC at high R. The tendency for the LSC to be in a preferred 
orientation due to the Coriolis force could be cancelled by a slight tilt of the apparatus relative to 
gravity, although this tilt affected other aspects of the LSC that the Coriolis force did not. 



I. INTRODUCTION 

The problem of Rayleigh-Benard convection (RBC) 
consists of a fluid sample heated from below 0, H 0|. 
In our case the sample was a cylindrical container filled 
with water. This system is defined by three parameters: 
the Rayleigh number R = agATL 3 /nv (a is the isobaric 
thermal expansion coefficient, g the acceleration of grav- 
ity, AT the applied temperature difference, L the height 
of the sample, k the thermal diffusivity, and v the kine- 
matic viscosity), the Prandtl number a = v/n, and the 
aspect ratio r = DJL (D is the diameter of the sam- 
ple). Convection of heat occurs primarily as a result of 
the emission of volumes of hot fluid known as "plumes" 
from a bottom thermal boundary layer that rise due to 
a buoyant force, while cold plumes emitted from a top 
boundary layer sink. In the turbulent regime of Y = 1 
samples that we studied, theseplumes drive a large-scale 
circulation (LSC) 0, HIE II 1111111113 which is oriented 
nearly vertically with up-flow and down-flow on opposite 
sides of the sample. 

While the LSC configuration does not have the cylin- 
drical symmetry of the sample, we found that the struc- 
ture underwent an erratic azimuthal motion that allowed 
it to sample all angles over time. We showed elsewhere 
that, on a long time scale (days), this azimuthal mo- 
tion could be partly described as a diffusive process. We 
believe it to be driven by the intense turbulent fluc- 
tuations that prevail throughout the system. Interest- 
ingly, relatively sudden re-orientations through various 
angles (with durations of order minutes) were superim- 
posed onto the diffusive motion at intervals much longer 
than their duration. These rc-oricntations could occur 
via two distinct mechanisms. One of these (a "rotation" ) 
consisted of a relatively fast rotation of the circulation 
plane (or of the vorticity vector of the LSC) without sig- 



nificant change of the circulation speed (or of the magni- 
tude of the vorticity vector) . The other (a "cessation" ) 
involved the complete decay of the circulation, followed 
by a spontaneous new start with a different orientation. 
These processes and their statistical properties have been 
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discussed in detail in previous work 

Over long time scales (days) we found that the LSC 
did not sample all angles equally, suggesting that the 
azimuthal diffusion was subject to some asymmetric po- 
tential. In addition, the meanderings of the LSC ori- 
entation occasionally added up to rotations through ap- 
proximately 27T, which we refer to as "revolutions" . For 
R > 10 10 clockwise (when viewed from above) revo- 
lutions were found to be more frequent than counter- 
clockwise ones. In this paper we report on these mea- 
surements, and develop a simple model for the effect of 
the Coriolis force on the LSC that explains many of these 
observations. Experimentally, the strongest qualitative 
evidence for the influence of the Coriolis force on the LSC 
is a change in the sample frame of the azimuthal loca- 
tion of the preferred orientation that occurred when the 
entire apparatus was rotated through a constant angle 7 
in the laboratory frame. 

In order to develop the model, we show that the 
Earth's rotation applies a torque to the horizontal flow 
components near the top and bottom plate that by itself 
and in the northern hemisphere would induce a rotation 
in the clockwise direction. In the absence of dissipation 
the torque in turn would be diminished by this rotation, 
leading to a rotation frequency precisely equal to that 
required for the torque to vanish and for the orientation 
of the LSC plane to be steady in an inertial frame. In 
the physical system one also needs to consider the dis- 
sipation experienced near the side wall and the top and 
bottom plates; we introduced this dissipation into the 
model by considering the viscous drag across a laminar 
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boundary layer and found that it reduces the rotation 
frequency. Next we need to introduce the coupling of the 
Earth's rotation to the vertical components of the flow 
near the side wall. It imposes a preferred orientation that 
is just north of west and that provides an azimuthally 
varying potential. This potential interferes with the ro- 
tation, and the two effects mentioned so far would impose 
a unique LSC orientation near the minimum of the po- 
tential. However, the turbulent background fluctuations 
that we have neglected up to now couple to the LSC and 
cause a meandering in the neighborhood of the poten- 
tial minimum, and induce occasional transitions across 
the potential barrier that result each time in a rotation 
through approximately one revolution in either direction. 
The clockwise and counter-clockwise transitions are not 
equivalent because the coupling to the horizontal veloc- 
ity components provides a slanting background to the po- 
tential (a "washboard" potential), leading to a preference 
for rotations in the clockwise direction when viewed from 
above, as seen in the experiment. Earlier [12J we showed 
that the LSC meandering has some qualities of a diffusive 
process, and derived its azimuthal diffusivity from mea- 
surements of the long-time angular fluctuations of the 
LSC orientations. We now use these diffusivities, which 
vary with R, in a stochastic model based on the potential 
discussed above to derive expressions for the probabil- 
ity distribution of the LSC orientation, for total number 
of rotations, and for the excess number of rotations in 
the clockwise direction from a Fokker-Planck equation. 
These results reproduce many of our experimental mea- 
surements. However, the model does not explain why the 
net rotation is absent for R < 10 10 . At large R we also 
found that, in addition to the Coriolis force, we had to 
invoke an additional azimuthal asymmetry (associated 
with small imperfections of the apparatus) in order to 
explain all of our data over wide parameter ranges. 

It is interesting to compare our results with experi- 
ments by Hart et al. where turbulent RBC samples 
of size comparable to ours were rotated deliberately at 
a constant rate in the laboratory frame about a vertical 
axis. That work involved rotation rates f2 much larger 
than the Earth's rotation, and therefore imposed much 
larger Coriolis forces upon the LSC. It showed that the 
LSC underwent an azimuthal precession in the frame of 
the rotating sample and in the retrograde direction at 
a rate which was proportional to Q when Q was not too 
large. We also find a net rotation of the LSC in the clock- 
wise (i.e. retrograde) direction. However, when scaled by 
the axial component of f2 its rate is larger by a factor of 
three or so. In this comparison it should be considered 
that the Earth's rotation rate of 7.3 x 10 -5 rad/s is more 
than two orders of magnitude smaller than the smallest 
experimental values used in Ref. |14|. A possible ex- 
planation for the quantitative difference may be found in 
a greater role of viscous dissipation at the larger fl values 
used in Ref. 0. 

Hart et al. also developed a model for the coupling 
of the Coriolis force created by their deliberate rotation 



to their LSC that in some ways is similar to ours. How- 
ever, in their case the coupling is only to the horizontal 
velocity component because the rotation axis is parallel 
to the vertical flow. Their model thus has no potential 
extrema, the rotation of the LSC is only in the retrograde 
direction, and stochastic effects need not be invoked to 
yield rotation. When we set the latitude equal to it/2 in 
our model and neglect stochastic effects, our model be- 
comes equivalent to that of Hart et al., except that those 
authors used a different approach toward modeling the 
dissipation. 

While the Coriolis force plays a major part in many 
geophysical and astrophysical convection systems, it is 
rare that it can be observed in laboratory fluid-flow ex- 
periments because of the small magnitude of the force on 
laboratory scales. One example is high-precision mea- 
surements of the velocity profile of laminar flow in a long 
pipe |15| . It was found that the expected parabolic profile 
of the velocity field was skewed by the Earth's rotation. 
In that case it was possible to carry out a quantitative 
calculation of this effect based on the equation of motion 
(the Navicr-Stokes equation) of the system because the 
unperturbed velocity was known analytically, the velocity 
field remained time independent, and the Coriolis- force 
perturbation could be calculated at linear order. For our 
system the turbulent flow field is much more complicated 
and a full treatment based on the equation of motion 
does not seem feasible except perhaps by direct numeri- 
cal simulation. Thus we resorted to the simplified model 
described above. 

A careful study of the effects of the Earth's Coriolis 
force could prove useful for understanding and analyzing 
data in other high-precision RBC experiments. For in- 
stance, our study made it possible to explore other asym- 
metries that affected the LSC in our experiment, in par- 
ticular effects due to our top-plate cooling-system and to 
small deliberate misalignments relative to gravity |l6j . 

In the next section we describe briefly the experimen- 
tal measurements; a more detailed account of the appa- 
ratus had been provided before 0]. Section fill Al gives 
the results accumulated over many months for the az- 
imuthal rotation of the LSC. Data for the diffusivity of 
the LSC meandering are given in Sect. IIII Bl The mea- 
surements of the probability distribution of the LSC ori- 
entation are discussed in Sect. IIII CI Measurements of 
the instantaneous rotation rate turn out to be useful to 
determine a parameter of the model, and thus are given 
in Sect. If f f Dl In Sect. II VI we present a model for the 
interaction between the Earth's Coriolis force, the LSC, 
and the background turbulence of the system. In a num- 
ber of subsections we derive predictions for several prop- 
erties of the LSC and compare them with experimental 
results. The properties include the preferred orientation, 
the effect of dissipation, the probability distribution of 
the LSC orientation, the net rotation rate, and the to- 
tal rate of occurrence of revolutions. Comparison with 
other experiments is provided as well. The influence of 
azimuthal imperfections in the sample cell, we believe 



3 



primarily due to the cooling system of the top plate, is 
incorporated in the model in Sect. [V] Section con- 
siders the influence of a slight misalignment relative to 
gravity of the sample, and discusses how this tilt can in 
part cancel the Coriolis-force effects. A summary and 
some conclusions are presented in Sect. IVHI 

II. EXPERIMENTAL METHOD 

The experiments were done with two cylindrical sam- 
ples with aspect ratio T « 1 that were the medium and 
large sample described in detail elsewhere 01 • Both had 
circular copper top and bottom plates with a plexiglas 
side wall that fit into a groove in each plate. The plate 
surfaces that touched the fluid sample were nominally 
symmetric and fine machined to a roughness of less that 
2 micro-m. The side walls had variations of their diame- 
ter of less than 10 -2 cm. There was a small semi-circular 
hole 1.6 mm in diameter at the edge of each plate that 
was used for filling and to allow expansion of the fluid. 
There were no internal sensors or other structures that 
could lock the flow direction. A resistive heater wire was 
embedded in the bottom plate in such a way that it cov- 
ered all regions of the plate equally for even heating. The 
top plate contained a double-spiral water-cooled channel 
which covered all regions of the plate, with inlet and out- 
let on opposite sides of the plate. The medium sample 
had dimensions D — 24.81 cm and L = 24.76 cm. The 
large sample had D — 49.67 cm and two side walls which 
could be interchanged with L = 50.61 and L = 49.54 
cm. Each apparatus was filled with de-ionized de-gassed 
water and the average temperature of the bottom and 
top plates was kept at 40.0° C where a = 4.38. The 
two samples of different physical size allowed us to cover 
a larger range of R at the same a and T, so the overall 
range studied was 3 x 10 s J; R ^ 10 11 . Both samples were 
carefully levelled to better than 0.001 rad, except for the 
experiments in which we deliberately tilted the samples. 

Eight blind holes, equally spaced azimuthally in the 
horizontal mid-plane of the samples, were drilled from 
the outside into each side wall. Thermistors were placed 
into them so as to be within d = 0.07 cm of the fluid 
surface. They were numbered i = 0, . . . , 7 in the counter- 
clockwise direction when viewed from above, with ther- 
mistor normally located in an easterly position except 
at times when the entire apparatus was rotated in the 
clockwise direction through an additional angle 7. Since 
the LSC carried warm (cold) fluid from the bottom (top) 
plate up (down) the side wall, these thermistors detected 
the location of the upflow (downflow) of the LSC by in- 
dicating a relatively high (low) temperature. No parts 
of the thermistors extended into the sample where they 
might have perturbed the flow structure of the fluid. 

We made measurements of the thermistor tempera- 
tures at the horizontal midplane at time intervals St ss 6 
seconds, and fit the function 

T, = T + 5coa(iir/4 - ),i = 0, . . . ,7 , (1) 
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FIG. 1: An illustration of a revolution. The figure shows a 
time series of the orientation 9o(t) for the large sample over a 
time interval of about 12 hours (1/2 day) at R = 8.9 x 10 10 . 

separately at each time step, to these eight readings. Ex- 
amples of such fits have been shown previously [T3 . Il3| . 
The fit parameter 5 is a measure of the amplitude of 
the LSC and 0$ is the azimuthal orientation of the plane 
of the LSC circulation in the sample frame. As defined 
here, the orientation 8q is on the side of the sample where 
the LSC is up-flowing and is measured relative to the 
location of thermometer zero. Typically the uncertain- 
ties for a single measurement were about 13% for 5 and 
0.02 rev. for 8 Q . When we fit the side-wall thermistor- 
temperatures to get 9 0l the fit only determined the value 
of #0 modulo 2ir. In order to track #0 continuously 
through many azimuthal rotations we chose the value 
that was within ir of the value measured at the previous 
time step. 

III. EXPERIMENTAL RESULTS 
A. Revolutions and net rotation of the LSC 

Elsewhere we presented measurements of sudden re- 
orientations of the plane of circulation of the LSC 
through a continuum of possible angles A9 0|. They 
were relatively rare events in that their collective du- 
ration took up only a small fraction of the observation 
time. They assumed the form of either cessations, where 
the flow actually stopped before re-starting in a new 
randomly chosen orientation, or of rotations where the 
flow continued but its vorticity vector rotated relatively 
rapidly through some angle A8. Rotations and Cessa- 
tions typically occurred rather quickly, in the time span 
of about one turnover time of the LSC (which is on the 
order of 100 s). Here we focus on yet another event which 
was like a rotation in the sense that the LSC orientation 
changed without a significant change of the LSC speed, 
but the change in orientation A9 was close to 2ir and re- 
quired a larger time interval, typically of an hour or so (of 
order 30 turnover times) . These events will be referred to 
as "revolutions" . An illustration of a revolution is given 
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FIG. 2: A time series of the orientation 6o{t) for the large 
sample over 18 days at R — 8.9 x 10 10 . There was a net 
rotation of 17 clockwise revolutions. 



in fig. n After a revolution the orientation (modulo 2ir) 
is essentially the same as it was before that event. In 
our previous analysis |l2j | revolutions were too slow to be 
counted as rotations. 

We shall show below that revolutions are associated 
with diffusive motion in a periodic potential with a min- 
imum at a particular azimuthal orientation. This po- 
tential is created to a large extent by the coupling of 
the Earth's Coriolis force to the LSC. It turns out that 
clockwise revolutions are more frequent than counter- 
clockwise ones, leading to a net rotation over long time 
periods. This asymmetry again can be explained in terms 
of the Coriolis force. 

An example of net rotation of the LSC for the large 
sample at R = 8.9 x 10 10 is seen in a time series of 0o(i) 
shown in Fig. [3 There is a net rotation of 17 clockwise 
revolutions over 18 days. Inspection of the time series 
shows that the LSC had a tendency to remain aligned 
near a preferred orientation m much of the time, but 
that it occasionally underwent significant positive or neg- 
ative rotations before locking in again near m (modulo 
2ir) as illustrated in Fig. ^ 

To facilitate automated data analysis, we define a pos- 
itive revolution as an event where 9o(t) rotates from 
the range 9 m — ir/4 < 9q < 9 m + ir/4 to the range 
(9 m + 2tt) - tt/4 < O < (9 m + 2tt) + tt/4. A nega- 
tive revolution is the analogous change of approximately 
— 2n. We require the orientation to come within a small 
range of 9 m to make sure that each revolution counted 
is a rotation of approximately 2tt and not jitter around 
an angle 9 m ± it. Let us describe the sequence of N rev- 
olutions by a sequence of N numbers r, which are equal 
to 1 (-1) for a counter-clockwise (clockwise) revolution. 
The net rotation 1Z is then equal to the sum of the N 
values of 7"j. For every adjacent pair of values r*j and rj+i 
we compute the product di — r^i+i. The parameter di 
will be equal to 1 if the two members of the pair of 
correspond to rotations in the same direction, and equal 
to -1 otherwise. The sum V of all N — 1 di is a measure 
of the correlation between successive events. 

Counting the number of revolutions over all of the ex- 



perimental running time in the large sample of 258 days, 
covering the range 3 x 10 9 5; R S; 10 11 , we found a net 
rotation 1Z — —77, corresponding to 77 clockwise revo- 
lutions out of a total of N = 439 revolutions. The net 
correlation between successive events over this time was 
T> = 1; the smallness of this number relative to the to- 
tal number of revolutions N shows that the directions 
of successive revolutions were not correlated with each 
other, and that the net rotation was not due to a long 
time scale over which the LSC rotated in one direction. 
It also allows us to use binomial statistics to analyze the 
significance of the results. If there were a 50% chance 
for each revolution to be either in the positive or neg- 
ative direction, we would expect that the net rotation 
after N = 439 revolutions would be less than about one 
standard deviation er/v = V^V which is approximately 21 
revolutions. Thus the observations of 77 clockwise rev- 
olutions, corresponding to 3.7 cr/v away from a mean of 
zero, indicate that the net rotation is a real effect with 
a confidence limit of 99.98%. The average rate of net 
rotation in the large sample was uj-ji = —0.30 ± 0.08 rev. 
per day (clockwise is negative in our coordinates) . In the 
medium sample, covering the range 3 x 10 s % R S; 10 10 , 
there were only 9 net counter-clockwise revolutions out 
of a total of 331 revolutions over 194 days, and T> = —8 
for an average rate of w-r = 0.05 ±0.09 rev. per day. One 
sees that the results from the medium sample are consis- 
tent with no net rotation, and at the least they indicate a 
smaller rotation rate than was found in the large sample. 

We compare our data to an experiment by Hart et al. 
fl4| with r = 1 and Prandtl number a — 8.4 that stud- 
ied the effect on the LSC of the Coriolis force due to 
rotating a convection apparatus at various applied rota- 
tion rates f2 over the range 0.015 rad/s % 57 ^ 0.7 rad/s. 
Two samples were used; one was about the same size as 
our large and another a little smaller than our medium 
sample. For their large sample the authors reported no 
net rotation of the LSC for R < 10 10 but a net rota- 
tion of lu-jz = — 0.23J7 in the rotating sample frame for 
i?~2.9-10 11 18]. This is seemingly consistent with 
our rotation rate of uj-jz = (—0.30 ± 0.08)0, even though 
their applied rotation rates were several orders of mag- 
nitude larger than our (i.e. the Earth's) rotation rate 
Q = 7.3 x 10~ 5 rad/s, although later we will show that 
such a straightforward comparison is probably not cor- 
rect. It is also interesting that neither experiment found 
net rotation in the smaller sized apparatus, which yielded 
data only for the range R $ 10 10 . 



B. Azimuthal diffusion of the LSC 

The re-orientations (cessations and rotations) were rel- 
atively rare and their duration occupied a relatively short 
fraction of a given long time series. During the remainder 
of the time the azimuthal orientation 0o meandered in a 
seemingly random manner. Elsewhere we reported that 
over long time scales the orientation meanderings have 
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FIG. 3: Time series of the orientation 8o(t) at R — 4.7 x 
10 10 for the large sample. The solid line is the original data. 
The dotted line is the result of removing a re-orientation by 
shifting the data beyond it along the horizontal and vertical 
axes. 



the diffusive quality that the root-mean-square (rms) ro- 
tation rate computed over various time intervals nSt is 
given by 



10° 10 1 10 2 10 3 
nSt (s) 




(2) 



FIG. 4: (a): The rms rotation rate #™ 3 as a function of 
nSt on logarithmic scales for the large sample and for R = 
5 x 10 9 . Solid symbols: from the analysis of the complete 
data sequence. Open symbols: from the analysis of the data 
sequence after removal of re-orientations. The solid line is 
a fit of Eq. |5|to the open symbols, adjusting De. (b): The 
Dg = 4.45 x 1CT 4 .R - 556 x v/L 2 jT^. Since re-orientations diffusivities D e , scaled by L 2 /u, as a function of R for the 



Here 9 r n ms = y/([0 o (t + nSt) - 9 {t)} 2 ) /(nSt), Dg is a dif- 
fusivity, St is the time interval between successive data 
points, n is a positive integer, and () indicates a time 
average We found an effective diffusivity given by 



do not occur very frequently, that analysis had used time 
series of 9o that included all re-orientations even though 
these events could not be described by simple diffusion. 
Motivated by the present need for more accurate val- 
ues of Dg, we repeated the analysis after filtering our 
re-orientations. The procedure is illustrated in Fig. 
There the solid line represents a short section of the orig- 
inal time series. It clearly reveals a re-orientation. The 
modified time series is shown as the dotted line. The dis- 
placements along both the time and the 9 axis during 
the re-orientation have been removed by shifts parallel to 
the axes. 

Figure 0^ shows an example of 9™ s as a function of 
nSt on logarithmic scales for both methods of analysis. 
One sees that omitting the re-orientations still yields data 
consistent with Eq. (solid line in the figure) , but that 
the diffusivity Dg is about a factor of two smaller when 
the re-orientations are omitted. The new analysis yielded 
the results for Dg shown in Fig. ^p. They can be repre- 
sented well by 



Dg = 2.3 x lQ-^R 



3 D 0.46±0.03 



x v/L 2 



(3) 



as is shown by the solid line. It is this diffusivity based 
on the modified time series for 9q that we use in the 
calculations presented in this paper. Presumably it bet- 
ter represents diffusive processes of 9o(t), as the values 



medium (solid circles) and the large (open circles) sample as 
determined after removal of re-orientations. 



reported earlier were " contaminated" by non-diffusive re- 
orientations. 



C. Preferred orientation of the LSC 

Over the course of this work, the entire apparatus was 
rotated through several different angles 7. We measured 
7 as a clockwise angle when viewed from above in the 
laboratory frame, with 7 = when 9 = was pointing 
east. After each such reorientation the top plate of the 
sample was carefully levelled to within 0.001 rad. A pair 
of time series 9o(t) before and after one such rotation is 
shown in Fig.JSJfor R = 9 x 10 9 in the large sample. The 
solid line corresponds to 7 = and the dotted one is for 
7 = 0.977T. While there were many fluctuations in 9o(t) 
for both time series, one sees that for each run there was a 
preferred orientation, and that this orientation shifted by 
about 7T in the sample frame when the sample was rotated 
by 7 ~ 7r, suggesting that this preferred orientation was 
due to an external field fixed in the laboratory frame and 
not due to an asymmetry fixed in the sample frame. 
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FIG. 5: Time series of the orientation 9o(t) for the large sam- 
ple at R — 9.4 x 10 9 for sample orientations 7 = (thick solid 
line) and 7 = 0.97-7T (thin dotted line). While there were many 
fluctuations of 0o(t), notice that for each run there was a pre- 
ferred orientation that the LSC tended to stay near much of 
the time, and that this orientation shifted by about n in the 
sample frame when the sample was rotated by n , indicating 
that the preferred orientation was due to an external field. 



We determined the preferred orientations 9 m of the 
LSC from the probability distributions of the orientations 
p(6*o). The Gaussian function 



P(0o) 



C_ 
2^ 



1-C 



27T<7g 



exp • 



1 1 



2ai 



(4) 



was fitted to the data for p{9q). Examples of such fits 
are shown in Fig. [St for the data of Fig. [5] and in Fig. 03 
for data with the same 7 = and 7 = 0.977T but with 
R = 9.4 x 10 10 . Although the coordinate 8 is periodic 
and the fitting function is not, this did not cause any 
difficulties since we found that the width of the Gaussian 
generally was much less than it. The only correction 
needed was a shift of the values of 9 by multiplies of 2tt 
so that 9q was within ±7t of 9 m . 

In Fig.[Hli the shift of 9 m [the peak of p(6q)} after ro- 
tating the apparatus is unmistakable. One can see that 
it is only a little smaller than 7. However, at larger R, 
the shift of p(9q) was much smaller than 7. This is shown 
in Fig. for R = 9.4 x 10 10 , for the same 7 = 0.97tt. 
The smaller shift suggests that the preferred orientation 
at this R is chosen in part by some sample asymmetry 
rather than exclusively by an external field. This issue 
will be discussed in detail in Sect. fVl 

The results for 6 m are shown in Fig.0 in both the sam- 
ple reference frame (a and c) and the laboratory reference 
frame (b and d). One sees clearly for the large sample 
in the sample frame (Fig. [7J:) that the data for differ- 
ent angles 7 approach each other at large R, indicating 
that there is an asymmetry associated with the appara- 
tus that dominates 9 m at large R. Comparing with the 
vertical dotted line in Fig. H3 one sees that the preferred 
angle in the sample frame at large R is close to the hot 
side of the top plate. This will be discussed in detail in 
Sect. E| 
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FIG. 6: (a): The probability distribution of the orientation 
p(6o) for the data shown in Fig. [£] The function p(0o) = 
C + (1 - C)/{V2^a e ) x exp[- (0 O - 6 m f/{2aj)} was fit to the 
data to get the preferred orientation m . Solid circles: 7 = 0, 
fitted by solid line. Open circles: 7 = 0.977T, fitted by dashed 
line, (b): p(0 o ) as in (a), but for R = 9.4 x 10 10 . The vertical 
dotted line corresponds to the "hot" side of the top plate (see 
section Sect. [3 for an explanation). 



For the large sample in the laboratory frame (Fig.[3l), 
we see that the data for different 7 nearly collapse at 
small R, indicating that there is an asymmetry associated 
with an external field, fixed in the laboratory frame, that 
causes 9 m to align to the west. Similar trends can be 
seen in the medium sample, but the data do not collapse 
well, suggesting that the asymmetries are weaker or more 
equally balanced there. 



D. The instantaneous rotation rate of the LSC 

The model for the influence of the Coriolis force on 
the LSC to be presented in Sect. IIVI involves one free 
parameter. Within the context of that model this pa- 
rameter can be determined from independent measure- 
ments of the time averaged instantaneous rotation rate 
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FIG. 7: The preferred orientation 6 m of the LSC vs. R for var- 
ious rotations in the laboratory frame of the samples through 
angles 7. (a): medium sample, sample reference frame, (b): 
medium sample, laboratory reference frame, (c): large sam- 
ple, sample reference frame, (d): large sample, laboratory ref- 
erence frame. 9 m = is east in the laboratory reference frame 
when 7 = 0. Solid circles: 7 = 0. Open circles: 7 = 0.257T. 
Solid triangles: 7 = 0.507T. Open triangles: 7 = 0.757T. Solid 
squares: 7 = 0.977T. Dashed line: west. Dotted line: the 
"hot" side of the top plate (see Sect. [V] below). 
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FIG. 8: (a): The average of the instantaneous azimuthal ro- 
tation rate (#o,i(#o)} as a function of 8q at R = 5 x 10 9 in the 
large sample. Solid line: a fit of (0o,i(0o)) — Asm 8q — B to the 
data, (b): The fit amplitude A vs. R. Solid circles: medium 
sample. Open circles: large sample. Solid line: model value 
A = 2a{lcos<j)/(12RI 1/2 + 1) with a = 5.0. 



(0 o ,i(0 o )) = ([0o(t + St) - 6 (t)]/5t) where () denotes an 
average over time but is left as a function of 9q. The 
purpose of the time average is to remove stochastic ef- 
fects so what remains represents a deterministic rota- 
tion rate. In Fig. [SJl we present experimental results for 
(0o,i(0o))- At each 8o it shows an average over many 
data points of 0o,i(0o) for R = 5 x 10 9 in the large sam- 
ple, where m was found to be fixed in the laboratory 
frame (see Fig. 01) and presumably determined primar- 
ily by the Coriolis force (as discussed in Ref. [l^. 6$^ 
was calculated with a slight time delay r after 0o was 
reached which satisfied the condition 8q — 0). Below in 
Sect. lIVBl we show that our model predicts the relation- 
ship (0o,i(0o)) = ^4 sin 0o + B. This function was fit to 
the data for (0o.i(0o)) by adjusting A [B is a known con- 
stant, see Eq.^]. The fit amplitude A vs. R is shown in 
Fig.|Sp, along with a solid line representing a model value 
for A derived from Eg. 1121 and to be discussed below. For 
larger R where the Coriolis force does not determine the 
preferred orientation, we inserted a second free param- 
eter into the fitting function to allow for a shift in m . 
This issue will be discussed in Sect. 
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IV. A MODEL FOR THE INFLUENCE OF THE 
EARTH'S CORIOLIS FORCE ON THE LSC 

A. Preferred orientation 

The Earth's Coriolis force can be shown to be responsi- 
ble for the net rotation of the LSC, and to dominate the 
preferred orientation for some parameter ranges. The 
contribution of the Coriolis force to the fluid acceleration 
in the sample frame is given by 



da 
~dl 




d 2 e Q 
dt 2 



(5) 



Here u is the fluid velocity, ft is the Earth's rotation 
rate of one revolution per day, and 9 = d0 /dt is the 
azimuthal rotation rate of the LSC in the frame of the 
sample. The sum ft + dOo/dt is the total rotation rate 
in an inertial frame. The centrifugal force is ignored be- 
cause it is of 2nd order in ft, which is a small parameter. 
We shall take cf> — 34.46°, the angle of latitude at our 
location, to be positive in the northern hemisphere. It 
is natural to consider u — (u r ,ue,u z ) in the cylindrical 
coordinates of the sample. 

Although the experimental data for the LSC orienta- 
tion 6>o in Figs. to [3] clearly show, on three different 
time scales, that this angle varies erratically in time, we 
shall ignore the stochastic effects for now as we derive the 
deterministic effects of the Coriolis force. If we assume 
that the deterministic forces are independent of the time- 
dependent stochastic forces, then the deterministic equa- 
tions are equivalent to time-averaged equations. Start- 
ing in Sect . IIV Dl below we will consider the influence of 
stochastic effects, due to the turbulent background fluc- 
tuations, on the flow. 

For simplicity we model the LSC as having a square 
vertical cross section and as consisting of four compo- 
nents: two horizontal for flow across the top and bottom 
plates and two vertical near the side wall. The flow is 
taken to be horizontal and flowing towards 9q in a cone 
based on the bottom plate and with its apex at the center 
of the cell. It is in the opposite direction in a congruent 
cone based at the top plate with its apex also at the cell 
center. The horizontal flow is assumed to have veloc- 
ity Uh — 2zU/L, where z is the axial coordinate mea- 
sured from the center of the cell, and U is the velocity 
at the peak of the profile near the thin boundary layer 
at z ~ ±L/2. The flow in the remainder of the sample is 
taken to be vertical with velocity u v — 2rU / L, and up- 
flowing at angles within ±7r/2 of 9q, and down- flowing 
at other angles 12] . The velocity scale U is much larger 
than r#0i so we can safely ignore the rotational contribu- 
tion to the flow velocities here. 

The horizontal components of the flow are deflected 
to the right in the northern hemisphere, so both the top 
and bottom legs of the LSC are pushed into a clockwise 
acceleration: 



We ignore forces in directions other than the azimuthal 
because they will be balanced by the wall, although they 
may contribute to some deformation of the LSC struc- 
ture. For the vertical components, the upfiow is deflected 
to the west while downflow is deflected to the east in the 
northern hemisphere. Since we defined the LSC orienta- 
tion so that there is upflow at 9q, the LSC will be pushed 
to orient 9o towards the west (9 = ir) according to 

ue,v{i~) = — 2u v (r)fl cos <fi sin(ir — 6*o) — 9^r 

= 2u v (r)S} cos 4> sin 9q — 9 r . (7) 

We can combine Eqs. El and [7| taking an average of 
the stated velocity profiles over the entire volume of the 
sample to obtain the average azimuthal acceleration. At 
this point we also introduce an adjustable parameter a 
to be used as a fitting coefficient, and obtain 

(u g )=U afl cos (f> sin 9 - (ft/2) sin - # /2 ~ s o(r) ■ 

(8) 

In doing this velocity average we assumed that u(r, 9, z) 
is not correlated with 9q or sin#o- The velocity average 
led to the second and third terms being half as strong as 
the first because the horizontal flow in the cone regions 
takes up 1/3 of the volume of the sample, and thus a 
volume half as large as that of the vertical flow region. 
The parameter a accounts for any other differences in the 
ratio of vertical to horizontal flow forces; if the model is 
at all reasonable, then a fit of its predictions to the data 
should yield a value of a that is of order one. By dividing 
by the spatially averaged radius (r) = L/3, and moving 
the other 9o to the left side, we can now write a single 
differential equation for 9q- 



(ug) 



4L 



2aft cos 4> sin 9q — ft sin <fi — 9q (9) 



ug,h(z) = -2u/j(z)(ftsin</> + 9 ) - 6» r 



(G) 



We resorted to doing spatial averages because the forcing 
is not linear in r and thus there will be some distortion 
of the LSC structure which we do not consider explicitly. 
We derived 9q as a characteristic angular acceleration 
representing the entire LSC. 

The first term in the bracket on the rhs of Eq. ori- 
ents the preferred orientation 9 m to the west, while the 
second term forces the LSC clockwise, moving the pre- 
ferred orientation to the north of west. In Fig. [7| the 
data for 9 m at various azimuthal sample orientations 7 
nearly collapse in the laboratory frame at small R in the 
large sample, and 9 m is found to be slightly north of 
west; so the model has successfully predicted the quad- 
rant of 9 m . Equation [5] has two equilibrium solutions, 
corresponding to 9 = 9 = 0, that yield 9 m . They sat- 
isfy sin9 m = (l/2et)tan0, of which the value of 9 m in 
the second quadrant (between north and west) satisfies 
d9o/d9 < and thus is a stable equilibrium. However, 
because of systematic uncertainties in the measurement 
of 9 m as well as because of possible influences from sam- 
ple imperfections that have not yet been included in the 
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model, we can not accurately obtain the value of a at this 
point. In the absence of noise we would expect the LSC 
to assume a uniquely defined orientation given by 9 m . 
We attribute the distribution of 9$ about 9 m shown by 
the data in Fig. to the action of the intense turbulent 
background fluctuations on the LSC. 

B. Effect of dissipation 

While Eq. [5] provides a differential equation for 9o, it 
does not take into account viscous drag, which we con- 
sider here. We make the assumption that the dynamics 
are strongly damped, so we can ignore the inertial term 
(9 = 0). Then we can predict the steady state azimuthal 
rotation rate 9q(9q) by assuming that the Coriolis force 
is balanced by viscous drag in a boundary layer of width 
A. The average deceleration (ug M ) of the LSC can be 
estimated as the drag in the boundary layer times the 
fractional volume of the boundary layer in the sample 
given by X/L for each plate and AX/L for the side for a 
total of 6X/L: 



(ug iV ) = (vV 2 u e ) w ^-j^Y 

« (10) 
L 

where in Eq. E3 w e expressed A ~ 0.5L(i? e ) -1 / 2 |2£| in 
terms of the Reynolds number R e = UL/v, To find the 
steady-state rotation rate 9q characteristic of the LSC, 
we equate Eq. [S] and Eq. ^| to obtain 

9 Q = Asm9 - B (11) 

where 



(shown in Fig [Ht-) to obtain the fit coefficient A = (4.8 ± 

0. 3) x 10 -4 rad/s. By equating this with the first term of 

Eci.lHlwe obtain a = Ax (12R~ 1/2 + D/(2ricos6) ~ 5.0. 
We see that a is somewhat larger than unity, but not 
excessively so given the many approximations made in 
the model. It is unclear which approximations would 
lead to such an underestimation of the preferred orien- 
tation forcing. The preferred orientation 9 m is given by 
s'm9 m = B/A = 0.98-7T and is consistent with those ob- 
tained from the peak of p(9q) shown in Fig.[7]for small R 
in the large sample. This confirms that the ^"dependent 
forcing is responsible for the non- uniform p(9$) in this 
parameter range. This agreement is significant because 
there was no free parameter in Ea. llll for determining 9 m 
(A is determined mainly by the sinusoid amplitude and 
has only a small effect on 9 m ). The fit amplitude A vs. 
R is shown in Fig. along with the model value for 
A from Eq.[Hwith a = 5.0, and with R e = 0.0345i? 1 / 2 
which was found in other independent measurements [l9| 
for the range of R studied here. The large scatter around 
the fit can be qualitatively attributed to asymmetries of 
the sample (see Sect. 0. The scatter also prevents us 
from accurately testing the i?-dependence of the model 
at this point. 

We note that the line in Fig. |SJd is nearly horizontal, 
indicating a very small R dependence of the model. This 
is because the denominator of A in Ea. ll2l is dominated by 

— 1/2 

the constant term equal to 1; the term 12R e ' ranges 
from 0.1 to 0.4 in the relevant range of R. Physically 
this implies that drag is relatively unimportant in slowing 
down the azimuthal rotation of the LSC, especially at 
higher R. Thus the rotation rate is limited mainly by 
the reduction of the Coriolis force in proportion to —6, 

1. e. by the rotation of the LSC itself in the retrograde 
direction. 



D. The probability distribution of the orientation 



2a£l cos0 
12Re 1/2 + 1 



and B 



D, sin < 



12R, 



-1/2 



1 



(12) 



As already seen above from Eq. |SJ the preferred orien- 
tation 9 m is found when 9q vanishes, given by sin# m = 
(I /2a) tan0 = B/A. 



C. Instantaneous rotation rate 

In order to determine the free parameter a of the 
model, we examine data for the time averaged instan- 
taneous rotation rate (#o,i($o)) given in Sect. HTTP] and 
equate it with 9q(9 ) of the model, which at this point 
still does not take into account stochastic effects and thus 
is equivalent to a time-average. At R = 5 x 10 9 in the 
large sample, where the Coriolis force was seen to deter- 
mine the preferred orientation, we fit Eq. 1111 to the data 



Calculating the net rotation rate Qn, the total num- 
ber of revolutions per unit time f2jv, and the probability 
distribution p(9q) from the model is more complicated 
because these properties are controlled by the turbulent 
"noise" that causes the LSC orientation to meander. In 
Sect. lIIIBl we reported that the orientation meanderings 
have the diffusive quality that the root-mean-square ro- 
tation rate 9™ s is given by 0™ s = y/Dg/(nSt), where 
Dq is a diffusion constant. We thus adopt the idea that 9q 
carries out a diffusive motion, driven by the background 
turbulence, in the potential created by the Coriolis force. 
From Eq. 1111 we find this potential to be given by 



V= 9 (9 a )d9 = A cos 9 a + B9 Q . (13) 



With no noise, the orientation would become locked in 
a local minimum of this washboard potential and there 
would be no net rotation. With strong noise, the local 



10 



maxima in the potential become insignificant and the net 
rotation rate is equal to — B. 

To calculate both p(0o) and the net rotation rate for 
an intermediate noise strength we follow the method of 
Ambegaokar and Halperin who solved a mathematically 
equivalent problem for the voltage across a Josephson 
junctions driven by a current source p]|. We start with 
a Fokker-Planck (or Smoluchowski) equation 



dp(0 o ) 
dt 



V[-p(e )6 (d ) + D e Vp(d )] . (14) 



The probability distribution evolves in time according to 
the sum of an advective and a diffusive probability cur- 
rent. This equation is valid in the case of strong damping, 
specifically where the diffusion velocity 3C/D/4L > 130 
is large compared to the mean drift velocity (less than 
rising), and the mean free path ALD/3U « 0.2 rad is 
small compared to the distance between peaks of the po- 
tential (2ir rad). In the stationary state the left side of 
Eg. 1141 is equal to zero, and by integrating the right side 
we obtain 



- p(0 o )0 o (9 o ) + D 6 Vp(0 o ) 



' 2ir 



(15) 



Here — w-r./(27t) is a constant of integration whose value 
can be confirmed by integrating Eq. 1151 once more over 
the range from to 2it and using the fact that p(0o) 
is normalized with periodic boundary conditions in this 
reduced range. In other words the net rotation rate is 
the balance of the advective and diffusive terms. 
Solving Eq. ^| for p{0) is done by assuming 
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FIG. 9: The probability distribution of the orientation p(9o) 
for R = 5 x 10 9 in the large sample. Open circles: data. 
Solid line: prediction based on the Fokker-Planck equation 
and the experimentally determined diffusivity De without any 
free parameters. 



large sample along with the result calculated from Eq.fPol 
using the experimentally determined values of Dg given 
in Fig. 21 in Sect. lIIIBl There is good agreement without 
requiring any adjustable parameters. Sample asymme- 
tries at larger R required the use of experimentally mea- 
sured values for A(R) in place of Eq . l*P2l and a shift in 0$. 
With these adjustments, similarly good agreement was 
found at all R studied in both samples. 



E. Net rotation rate 



P(0o) 



2nD e 



f(0o)g(0o) 



where 



f(0 o ) = exp{-V/De) 



(16) 



(17) 



Here p(0a) oc f(0a) happens to be the solution for B = 0, 
and we note that Eq. IT3l implies —VP = f?o- We apply 
periodic boundary conditions and get 



g(0 o ) 



/(27T) ' /(0) 



/(0) 



e ° dff 



f(O') 



/(27T) 



d9' 



ie /m 

(18) 

We next integrate Eq. Hoi from to 2ir. Requiring p(6o) 
to be normalized, we then have 221 



ujtz = 2irD g 



d0 o f(0 a )g(0 o ) 



(19) 



These equations must be solved numerically to obtain lo-ji 
and p{9q). Figure shows p(0o) for R — 5 x 10 9 in the 



Equation ^5] can be used to numerically calculate the 
net rotation rate uj-ji. The experimental data for are 
shown in Fig. 110b (circles). The predictions from the 
model with a = 5.0 in Eq. ^| and Dg from Eq. [3] are 
shown as well (lines). The agreement is less than per- 
fect. Since the Coriolis force was apparently not the only 
asymmetry contributing to the potential barrier V", we 
also calculated lo-k by substituting the measured values 
of A(R) and Dg(R)ior Eq. 1121 and obtained the triangles 
in Fig. I10fa . The agreement with the measurements (cir- 
cles) is improved somewhat. In either case, the model 
predicts lu-r close to — with some variation with R 
for both samples. While there is overall consistency with 
the experimental data from the large sample, reproduc- 
ing correctly both the direction and magnitude of the net 
rotation, the model does not explain why we see no sta- 
tistically significant net rotation in the medium sample. 

It is interesting to note that the prediction of \lu-ji\ in- 
creases with R. This does suggest that there is a lower 
cutoff for R below which net rotation cannot be observed, 
with the cutoff depending on the resolution of the exper- 
iment. For example at R — 5 x 10 9 in the large sample, 
the prediction is lu-ji = —0.03, which is too small for us to 
distinguish from zero. However, in our experiments cells 
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the model we follow the method of Kramers [2j|. He 
calculated the rate of crossing a potential barrier for large 
barriers. While in our case the diffusion constant Dg is of 
the same order as and sometimes larger than A, we can 
still calculate the ratio lun/lu-jz for these intermediately 
sized potential barriers. 

In our model 8q diffuses in a washboard potential with 
a series of minima. We would like to make use of the 
Fokker-Planck equation for the non-equilibrium proba- 
bility distribution p m (6o) of the LSC trapped in the m th 
potential well. While our system is not in the large- 
potential-barricr limit of the Arrhenius-Kramers prob- 
lem, the barriers are large enough to keep the LSC in a 
single well on average for several hours (2it/ujn) as can 
be seen, for example, in Fig.0 We claim that the barri- 
ers are large enough so that p m (&o) is nearly stationary 
and the stationary Fokker-Planck equation (Eq. I15|l is 
a good approximation in this case. Since p m (#o) is not 
periodic, we can treat two potential barriers AV+ and 
AV- (which have different sizes because of the net ro- 
tation term B in the potential) separately to obtain the 
rates ui± of transition from one potential minimum m to 
an adjacent minimum m±. We rearrange Eg. 1151 for a 
stationary distribution, replacing p(9o) with p m (#o) and 
Ufi with uj± to obtain 



FIG. 10: (a):The net rotation rate ujjz of the LSC vs. R. 
Solid circles: medium sample data. Open circles: large sam- 
ple data. Solid line: model calculation of uj-r. with A given 
by Ea. 1121 and a = 5.0 for the medium sample. Dashed line: 
model value of for the large sample. Solid triangles: cal- 
culation of uj-r, for the medium sample where the potential 
barrier A(R) is obtained experimentally at each R as dis- 
cussed in Sect. ITTTTjl Open triangles: calculation of uj-r. for 
large sample with experimentally obtained A(R). Because net 
rotations are infrequent, the data for uj-r (circles) were aver- 
aged over multiple runs in a small range of R. (b): The rate 
of revolutions ujn vs. R. Symbols have the same meaning as 
in part a. The calcuations of ujn are from Eq. 1231 



of different heights were used to cover different ranges 
of R, and because of the factor of L 2 in D, increasing 
R by increasing L (and keeping AT constant) leads to 
an increase in oj-ji with decreasing R in our parameter 
range. Thus even though the theory predicts a cutoff 
R below which net rotation cannot practically be ob- 
served, it does not explain the absence of net rotation in 
our medium sample. The limit of the net rotation rate 
as Dg/A becomes large at large R in both samples is 
= -B ~ -0.5ft. 



F. Frequency of revolutions 

Net rotations are very rare and thus the experimental 
uncertainty of their rate of occurrence is relatively large. 
A stricter test of the model is achieved by calculating the 
total number of revolutions ljm- To calculate ujn from 



f(0 o )uj ± « -2nD g V\p m (0o)f(O o )] 



(20) 



for #o between m and m±. Next we integrate over the 
potential barriers from m to m± to obtain 



UJ± 



2irD, 



M0o)/(0o)] 



mi 



(21) 



We note that f(9o) has the significance of the Boltzmann 
factor in the Kramers approach. We do not attempt to 
evaluate p m (#o) f° r small potential barriers because we 
only need to evaluate the ratio lu+/uj-. Since p m {0o) 
represents a particle trapped near m, it is negligible near 
m±. Now the numerator has the same magnitude in 
both the + and — cases and we can calculate the ratio of 
transition rates. Making use of Eas.lT7land[T51to evaluate 
the ratio of the remaining integrals we obtain 



UJ- 



C Wo + 2ir)d0 o 



— exp • 



-271-5 



(22) 



With this the ratio lun/lu-ji is found to be 



u + - uj_ D g 

We show the data for ujn in Fig. llOb in a similar man- 
ner to the data for lj-jz. The model calculation for u>n 
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(lines) uses the value of from Eq. ^5] with A obtained 
from Eq. The model captures the general trends of 
the data well. When the experimentally obtained A(R) 
are used in place of Eq. ^| (triangles) there is excellent 
agreement between the measured values of lon and the 
model predictions in both samples. 

Since the potential tilting parameter B increases the 
size of AV+ while equally decreasing the size of AV-, the 
value of B has a negligible effect on wjv- On the other 
hand the difference betwene AV+ and AV- is responsible 
for the net rotation, so B has a significant effect on ajfc. 
Since the model predicts the correct Un for both samples, 
the disagreement between the prediction of oj-ji and the 
data in the medium sample suggests that the values of 
B we used for the medium sample are suspect. As seen 
from Eq. the value of B is essentially independent 
of experimental measurements. Thus the disagreement 
for the medium sample does not imply a shortcoming 
of the Fokker-Planck description but is an indication of 
unexpected physics. 

We note that the Kramers approach relies on the exis- 
tence of (preferably large) potential barriers. According 
to our model, the potential barriers disappear completely 
when B > A, which is the case for the Earth's Coriolis 
force near the north or south pole, or for intentionally 
rotating experiments. In these cases our calculation of 
ljn would not apply. 



G. Comparison with other work 

Another experiment [24j similar to ours, at (f> = 22°, 
r = 1, a w 5, 10 9 < R < 10 10 , and L = 19.5 cm, 
found a clockwise net rotation at a rate that increased 
with R, which we estimate as an average over all R rep- 
resented in Fig. 4 of Ref. to be lj r « (-1.3 ± 0.4)J1. 
Although the runs were much shorter and thus the un- 
certainties greater than ours, the net rotation rate re- 
ported was larger than any found by us, even though 
the measurements were in the range of R and L where 
we found no significant net rotation. More alarming is 
that the magnitude of the rotation rate was larger than 
J7sin0 ~ 0.3751, the maximum allowed by the Coriolis 
force at this latitude. Since the time series of the ori- 
entations showed a strongly locked preferred orientation, 
one would expect that the potential barriers would have 
significantly reduced lu-jz from the azimuthal drift rate, 
as was the case in our experiment. We conclude that our 
model is not consistent with the data of Ref. j24[ , and it 
does not seem likely that any Coriolis-force-based model 
could account for the large rotation rate observed in that 
work. 

An additional similar experiment from the same labo- 
ratory |2H at R = 5.3 x 10 10 , a = 5.3 but with aspect 
ratio r = 0.5 was seen to have a much faster average net 
rotation rate of uj-jz = +13.551. This flow was visualized 
to show that there was a single LSC roll, so our model 
should apply to that experiment with some minor mod- 



ification. In this case the absolute rate of rotation was 
very large, and the rotation was counter-clockwise when 
viewed from above whereas a Coriolis-force model would 
predict a much smaller rate in the clockwise direction for 
any experiment in the northern hemisphere. These facts 
suggest that it is not the Earth's Coriolis force that is 
responsible for the net rotation in that case. 

We can compare the predictions of the model with 
the results of Hart et. al. 14]. They deliberately ro- 
tated their samples in the laboratory frame. Since their 
axis of rotation was along the central vertical axis of 
the cell and their rotation rate was much larger than 
one revolution per day, effectively <j) — tt/2 and the 
Earth's Coriolis force may be neglected. Since the ro- 
tation rate was relatively fast, one would expect any 
potential barriers to be less relevant, and in any case 
those due to the Coriolis force vanish for <fi = tt/2. 
Thus our calculation applicable to their system is that 
of the azimuthal drift rate with sin(</>) = 1 which gives 

lo n = uj n = -B = -n/(12Re 1/2 + l) ~ -0.951. The sign 
of this prediction agrees with but the magnitude is larger 
than that of the experimental result uj-jz = — 0.2311. Hart 
et. al. did not see any net rotation below R = 9 x 10 9 , 
which they achieved with a smaller sample about the 
same size as our medium sample. This is in striking sim- 
ilarity to our own observations, and suggests there is a 
physical reason related either to the size of the sample or 
to R that no net rotation was observed in these experi- 
ments. 

Hart et al. also developed a model to explain their 
results that is similar to ours in that the Coriolis force 
was balanced by drag. In the limit of zero drag, called the 
"inertial oscillation" model by Hart et al., both models 
yield lo-jz = —51. Better agreement with the experiment 
can be achieved for a dissipation level significantly larger 
than that provided by the viscous drag across a lami- 
nar boundary layer assumed by us. Indeed Hart et al. 
invoked an effective turbulent viscosity and introduced 
a greater level of dissipation into their model, thereby 
achieving good agreement with their measurements of 



V. EFFECT OF AZIMUTHAL SAMPLE 
ASYMMETRY ON THE PREFERRED 
ORIENTATION 

At large R we found a preferred orientation that was 
locked in the sample reference frame (see Fig.[7| and thus 
must have been due to some asymmetry of the apparatus. 
We identified the top-plate cooling-system as the most 
likely source. This system consisted of a channel in the 
form of a double spiral with inlet and outlet at opposite 
ends of a diameter. As coolant traversed through the 
channel, its temperature increased slightly as it absorbed 
heat from the plate. This meant that the coolant was 
coolest near the inlet and warmest near the outlet. We 
confirmed this by conducting separate experiments with 
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coolant flowing in opposite directions. In both cases the 
temperature monitored at five locations in the top plate 
near the fluid-solid interface showed that here was a small 
lateral temperature gradient in the top plate that was 
positive from the inlet to the outlet. In both cases the 
preferred orientation 9 m at large R corresponded to the 
warmer side of the plate. 

The double-spiral design had been intended to pro- 
vide a uniform temperature over the entire top plate be- 
cause ideally the average temperature of adjacent chan- 
nels should be the same everywhere. However, the ex- 
perimental results suggest that this cancellation between 
adjacent channels was not sufficient to prevent a biased 
azimuthal orientation of the very sensitive LSC. We be- 
lieve that the incomplete cancellation near the inlet and 
outlet may be responsible for the problem. The cooling 
of the top plate contrasts with the heating of the bottom 
plate, where the same amount of heat is dissipated per 
unit length along the entire length of the heater wire and 
where the spacing between adjacent lengths of the wire 
was smaller than that between the spiral channels. 

In the remainder of this section we model the effect 
on the LSC orientation of a small horizontal tempera- 
ture gradient in the top plate due to the cooling-system 
geometry. First we focus on the horizontal temperature 
different ST in the top plate between two thermistors 
a distance 3L/4 apart along a diameter nearly passing 
through the inlet and outlet. The magnitude in the large 
sample for R > 2 x 10 10 was given approximately by 
ST/ AT = 1.2 x Kr 8 i? a57 . A top-plate thermistor tem- 
perature and the side-wall thermistor temperature below 
it showed a positive correlation, with the top-plate tem- 
perature leading the side-wall temperature, suggesting 
that the top-plate temperature determined the LSC ori- 
entation and not the other way around. For smaller R the 
temperature gradient fell below the fit and even became 
slightly negative from the inlet to the outlet. At these R 
the Coriolis force was shown to determine 9 m , and the top 
plate happened to be aligned so that the cooling system 
pushes the LSC in the opposite direction as the Corio- 
lis force, so the temperature measurements suggest that 
the LSC determined the top-plate temperature-profile at 
small R. 

Next, to model the effect of the top-plate asymmetry 
on the LSC, we will assume that the top-plate temper- 
ature varies linearly along a diameter so that its mag- 
nitude is given by Ar5T/{3L)\ cos(6> - O )|- If this ST is 
added to the thermal boundary-layer temperatures, then 
it adds an additional buoyant force to the boundary layer 
of which a component proportional to sin(# c — 9q) will 
orient the LSC towards 9 C with the forcing 
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Now we take an azimuthal average of the acceleration 
over the entire sample by doing the integral of iig over 
9 to get a factor of 2/tt. Multiplying by the fractional 
thermal boundary-layer width 1/ (27V) one obtains 
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FIG. 11: (a): The preferred orientation 6 m vs. R in the 
laboratory frame for various rotation angles 7 of the large ap- 
paratus. Solid circles: 7 = 0, fit by solid line. Solid triangles: 
7 = 0.507T, fit by dashed line. Solid squares: 7 = 0.97-7T, fit by 
dotted line. The fitting function represents 9, n due to a combi- 
nation of the Coriolis force and the top plate asymmetry, (b) : 
Much of the same data plotted as 9 m vs. 7 for several values 
of R in the frame of the large sample. Points represent data, 
while the lines represent the model calculation using the fit 
values from Fig. a. Solid circles and solid line: R — 1 x 10 . 
Open squares and dashed line: R = 4.7 x 10 10 . Solid trian- 
gles and dotted line: R = 9.4 x 10 9 . Open diamonds and 
dashed-dotted line: R = 5 x 10 9 . 



the angular acceleration 



yo = 
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(24) 



Using the results TV = 0.0602i? 1 / 3 [HI and R e = 
0.0345-R 1 / 2 of prior measurements for our samples, 
we see that the acceleration scales roughly as R 12A , com- 
pared to the Coriolis- force acceleration from Eq.[5] which 
scales as R 5 ; so this explains why the Coriolis force de- 
termines 9 m at small R but the cooling system determines 
9 m at large R. To find the new preferred orientation, we 
set the forcings given by Eq. [5] and Eq. [2] equal to each 
other, setting 9q = 0, and inserting 7 in the Coriolis-force 
term to account for the different apparatus orientations. 
We obtain a relationship between R and 9 m given by 
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Here Rt is the transitional Rayleigh number where the 
cooling-system asymmetry overtakes the Coriolis force 
in determining 8 m . From the numbers given we calcu- 
late R t ~ 3.6 x 10 10 . We fitted our data for 9 m vs. 
R from the large sample for several apparatus orienta- 
tions 7 to Eq. [23 This fit is shown in Fig. II lb and gave 
R t — 3 x 10 10 and C = 0.267T. With these parameters, 
we also calculate 9 m vs. 7 from the model for several 
values of R. The results are shown in Fig. II lb for a 
different perspective on much of the same data as those 
shown in Fig.lllb. The transitional Rayleigh number R t 
is predicted quite accurately for an order-of-magnitudc 
model. The large-i? limit of the preferred orientation 
9 C is opposite the coldest top-plate thermistor, which is 
the first thermistor the fluid passes in the cooling chan- 
nel, so these data are consistent with the cooling sys- 
tem being responsible for determining 9 m at large R in 
our experiments. The cooling system orientation 6 C was 
shown in Fig.|U>. In that figure, the background of p(6q) 
for the 7 = data is notably larger than for the other 
p(6o) shown, which can be explained by the fact that at 
7 = 0, the Coriolis force and the top-plate asymmetry 
are trying to force the LSC into two opposing preferred 
orientations and the partial cancellation of these forces 
results in a more uniform p(9q). The variation of the po- 
tential barrier A with R as shown in Fig. |SJd can also be 
explained qualitatively. Since the cooling-system forcing 
and the Coriolis force oppose each other, the net effect 
is to reduce the net potential barrier at large R in the 
large sample from the Coriolis force prediction. A hor- 
izontal temperature gradient across the top plate from 
the medium sample was also measured at large R to go 
from the coolant outlet to the inlet. In this case the top 
plate was placed with the inlet on the east side, so that 
the cooling system forced the LSC in the same direction 
as the Coriolis force. This can be seen to increase the po- 
tential barrier A relative to the Coriolis force prediction 
for the medium sample in Fig. [SJd. 
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FIG. 12: The probability distribution of the orientation p(9o) 
for several tilt angles at R = 9.0 x 10 9 in the large sample. 
Dashed line: ,3 = 0. Solid line: (3 = 0.0035 rad. Dotted line: 
(3 = 0.0070 rad. Dashed-dotted line: (3 = 0.0105 rad. Notice 
that p(6o) is most nearly uniform for (3 = 0.0035 rad. 



parameter 


R = 9.0 x 10 9 


R = 8.9 x 10 10 




0.0026 ± 0.0001 


0.0021 ± 0.0001 




0.0031 ± 0.0002 


0.0014 ±0.0001 


A 


0.0029 ± 0.0004 


0.0015 ± 0.0004 


5 


-0.0044 ± 0.0009 


0.0056 ± 0.0009 


LO r 


-0.0032 ±0.0011 


0.0072 ± 0.0010 




-0.0031 ± 0.0021 


0.0058 ± 0.0039 


Re 


-0.0023 ± 0.0010 


0.0058 ± 0.0018 



TABLE I: The values of the tilt angle /3o where the probability 
distribution of various measured parameters (left column) had 
a maximum for the large sample. Results are given at two 
values of R. The preferred orientation at (3 = is determined 
by the Coriolis force at the smaller R, and by asymmetries of 
the top-plate cooling-system at the larger R. 



VI. CANCELLATION OF PREFERRED 
ORIENTATION BY TILT 

We can tilt the apparatus to an angle j3 relative to the 
gravitational acceleration, defined to be positive when 
the usually cooler side of the top plate is raised so as 
to add buoyancy to the boundary layers to oppose the 
prevailing preferred orientation. At an appropriate angle 
[3q > 0, we expect the effects of the various asymmetries 
to be cancelled out by the added buoyancy due to the 
tilt, ft would be a simple matter to combine the mod- 
els presented above for the Coriolis force and cooling- 
system force with the tilt model presented in Ref. [Ig 
to estimate /?o, but in the interest of brevity we show 
only the experimental results. The probability distribu- 
tion p(9o) is shown for several small positive tilt angles 
in Fig. El at R — 9.0 x 10 9 in the large sample. These 
experiments were done at 7 = 0, where the Coriolis force 



and cooling-system asymmetries oppose each other, and 
the Coriolis force dominates at the smaller R while the 
cooling system dominates at the larger R. For larger 
tilt angles p(0 ) is sharply peaked near the southeast 
(6 /2tt = 0.875), which corresponds to the high side of 
the tilted sample, as we would expect. At the interme- 
diate tilt angle of j3 = 0.0035, we find the most uni- 
form p(9o), suggesting a close balance between the initial 
asymmetries of the sample and the buoyancy effects due 
to tilting. To get a precise value of /3q, we fit Eq. 0]to 
p(0q) to obtain the width <jg and background constant C. 
Each fitparameter x was then fit to an empirical function 
x (j3) = xo+xx exp(-|/3-/3 |/x 2 ) for \(3\ < 0.021 rad. The 
values of (3q from the fit are shown in Table 1. 

Several other parameters were measured at various tilt 
angles to find (3q at two values of R in the large sample, 
with the results for {3 shown in Table 1. We previously 
reported the frequency of reorientations (spontaneous 
changes of the LSC orientation) oj r vs. (3, and showed 
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that a Gaussian function fit to u> r (0) with a maximum at 
some (3 Ql which has the same physical meaning as before. 
We found fa = 0.0093 ± 0.0010 rad at R = 9.4 x 10 10 
in the large sample 0| and (3 = 0.0022 ± 0.0006 rad 
at R = 1.1 x 10 10 in the medium sample 01 The 
other parameters measured were the coefficient A repre- 
senting the potential barrier strength; the LSC temper- 
ature amplitude 5; the root-mean-square rotation rate 
9rms', and the Reynolds number R e . In each of these 
cases the parameter x' is fit to an empirical function 
x'(J3) = x' + x[\f3 - /3 \ for \/3\ < 0.021 rad, with the 
value of Po again representing the balance between tilt 
and other asymmetric forces. Since ft is defined posi- 
tive when the tilt-induced buoyancy opposes the prevail- 
ing preferred orientation of the LSC, a positive value (3q 
means that the prevailing asymmetries affect the mea- 
sured parameter as expected. At R = 8.9 x 10 10 , all of 
the values of (3q are positive, so the cooling system asym- 
metry affects all of the measured parameters. This is 
reasonable because the tilt model and the cooling-system 
model both assumed the mechanism of increased buoy- 
ancy in the boundary layer. However, at R = 9.0 x 10 9 
some of the values of fio are negative, indicating that the 
prevailing asymmetry, in this case the Coriolis force, does 
not affect those measured parameters as strongly as the 
cooling system - which also opposes the Coriolis force 
and thus enhanced the tilt effect at this R - even though 
the Coriolis force had a stronger effect in determining 9 m 
at this R. The values of (3q are smaller in magnitude for 
the parameters relating to p(9o) at both R because the 
Coriolis force and cooling system asymmetries both apply 
and oppose each other, so there is a partial cancellation 
of forces. 

These data from the tilting experiments further con- 
firm that the Coriolis force does affect p(6q) and the local 
azimuthal forcing (<?o(#o)) that is responsible for p(0o), as 
we had already shown by other methods, but they also 
show that the Coriolis force does not significantly affect 8, 
R e , u) r , or 9 rms . Because the Coriolis force only deflects 
moving fluid and does not enhance or damp its speed, it is 
entirely reasonable that it does not directly affect R e and 
S. The other values, oj r and 9 rms , are presumably due to 
turbulence, but we do not know the detailed mechanisms. 
We can not rule out the possibility of weaker, "higher or- 
der" couplings between the Coriolis force and the LSC. 
At R — 8.9 x 10 10 , the cooling system increases R e by an 
amount equivalent to a tilt of /? = 0.006 rad. We previ- 
ously reported that R e increases by 185% per rad of tilt 
[l6| . so our results suggest that the cooling system has 
the effect of increasing R e by about 1% at this R. We 
measured Af as well, but our resolution of 0.1% is not 
enough to measure such small differences. We have pre- 
viously reported the sensitivity of Af to tilt is a reduction 
of only 3.1% per radian of tilt |16| . so even if we expected 
to find fio — 0.006 rad, that would imply a reduction in 
Af of only 0.02%. However, we consider even this small 
effect from the Coriolis force unlikely, as much evidence 
exists that Af is determined by the thermal boundary 



layers and not the LSC. 



VII. CONCLUDING REMARKS 

We identified several interesting effects of the Earth's 
Coriolis force and the turbulent background fluctuations 
on the large-scale circulation of turbulent Raleigh-Benard 
convection. There are clockwise and counter-clockwise 
relatively sudden revolutions through 2tt of the orienta- 
tion of the LSC circulation plane. The clockwise ones are 
slightly more frequent than the counter-clockwise ones, 
leading to a net rotation with an average rate of less 
than one rotation per day. There is a preferred value 
and a distribution about it, sampled over time, of the 
LSC orientation. Both the preferred value and the dis- 
tribution can be understood quantitatively in terms of 
the diffusive LSC meandering in the Coriolis-force po- 
tential, driven by the turbulent fluctuations. This model 
agrees with experimental results extremely well for small 
Rayleigh numbers in one of our samples, but at larger 
R a small asymmetry of the experimental apparatus has 
to be invoked to explain all the data. We gave a more 
complete summary of our findings in the Introduction 
and will not repeat this here. Instead, we make a few 
additional observations. 

We used two samples: a "medium" one over the 
Rayleigh-number range 3 x 10 8 ^ R ^ 10 10 and a "large" 
one over the range 5 x 10 9 S; R ^ 10 11 . We saw no net 
rotation in our medium sample, and one might wonder 
why this was the case. We suspect that the reason may 
be a change to a different flow structure with decreasing 
R that responds differently to the Coriolis force. A dif- 
ferent flow structure is sug gested based on images of the 
LSC reported by Ref. |2g. Their experiments are done 
with L = 1 and a = 4, with images taken by particle im- 
age velocimetry at R = 3.8 x 10 9 (which we can achieve 
in our medium sample) and R = 3.5 x 10 10 (which we 
can achieve in our large sample). The LSC at the higher 
R is more square in shape, filling out the container, while 
the LSC at the smaller R more nearly has the shape of 
an ellipse oriented with its long axis along a diagonal of a 
vertical sample cross section. In the latter case there are 
counter-rotating eddies at opposite corners of the sample. 
Since these smaller eddies rotate in the direction opposite 
the LSC, we would expect that the Coriolis force would 
have the opposite effect on them, which would partially 
cancel out the net effect on the LSC, and resulting in 
less net rotation than otherwise expected. These eddies 
would also be expected to reduce A and thus widen p(9 ) 
and increase lun, but since our model calculations used 
the experimentally obtained A(R) this does have any im- 
plications for the applicability of the Fokker-Planck de- 
scription. 

The calculation of p(9o) and the frequency of revolu- 
tions u>n from the Fokker-Planck equation, the Coriolis- 
force potential, and the turbulence-driven diffusion is re- 
markably accurate, which suggests that treating the LSC 
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meanderings as a diffusive process in a potential is an 
excellent model for the system. However, this does not 
completely describe the azimuthal dynamics of the LSC 
because it does not account for other spectacular events, 
in particular rotations and cessations. We found previ- 
ously by numerical simulation that a diffusive model un- 
derestimates the frequency of reorientations, even with- 
out potential barriers 01 • Cessations of the LSC have a 
duration of about one turnover time of the LSC, and rota- 
tions have a duration of less than 10 turnover times |12| . 
The potential barriers identified arc too weak to have a 
significant effect over such short time scales. The fact 
that Dg, from time series with reorientations removed, 
leads to an accurate calculation of ujn suggests the per- 
haps surprising result that rotations and cessations do 
not contribute significantly to lon- If we had used Dg 
with reorientations included, then we would have over- 
estimated ujn by about a factor of two. 

Studies of the asymmetry of our top-plate cooling- 
system provided some useful information for future ex- 
periments. While the LSC orientation and related dy- 
namics have proven to be very sensitive to asymmetries 
of the system in general, we have shown that several as- 
pects of the LSC, including R e , 5, 6 rms , and u> r , are not 
significantly affected by the Coriolis force, while they are 
affected by cooling-system asymmetries. In many stud- 



ies, the latter group of properties is a more physically 
interesting aspect of the LSC than p(8o), so while the 
Coriolis force may not affect the results of these stud- 
ies, a poorly designed cooling system will, particularly at 
large R where the temperature gradients in the cooling 
system become larger. The problem we had was a cool- 
ing system in which the temperature drop of the coolant 
as it flowed through the cooling channel left an associ- 
ated thermal imprint on the plate, which then affected 
the LSC temperature profile. To minimize this effect, the 
temperature drop in the coolant - given by the ratio of 
heat dissipated by the coolant to the total heat capacity 
of the coolant - should be small compared to the LSC 
temperature amplitude. Alternatively, the cooling sys- 
tem should be designed with such a geometry that the 
temperature drop of the coolant does not cause a hori- 
zontal temperature gradient in the plate. 



VIII. ACKNOWLEDGMENTS 

We are grateful to Werner Pesch for stimulating discus- 
sions. This work was supported by the US Department 
of Energy through Grant DE-FG02-03ER46080. 



[1] See, for instance, E.D. Siggia, "High Rayleigh number 
convection", Annu. Rev. Fluid Mech. 26, 137 (1994) for 
early work. 

[2] L. Kadanoff, "Turbulent heat flow: Structures and scal- 
ing", Physics Today 54, 34 (2001). 

[3] G. Ahlers, S. Grossmann, and D. Lohse, " Hochprazision 
im Kochtopf, Physik Journal 1, 31 (2002). 

[4] B. Castaing, G. Gunaratne, F. Heslot, L. Kadanoff, A. 
Libchaber, S. Thomae, X.-Z. Wu, S. Zaleski, and G. 
Zanetti, " Scaling of hard thermal turbulence in Rayliegh- 
Benard convection", J. Fluid Mech. 204, 1 (1989). 

[5] S. Ciliberto, S. Cioni, and C. Laroche, "Large-scale flow 
properties of turbulent thermal convection", Phys. Rev. 
£54, R5901 (1996). 

[6] D. Funfschilling and G. Ahlers, "Plume motion and large- 
scale circulation in a cylindrical Rayleigh-Benard cell", 
Phys. Rev. Lett 92, 194502 (2004). 

[7] R. Krishnamurty and L.N. Howard, "Large-scale flow 
generation in turbulent convection", Proc. Nat. Acad. 
Set. USA 78, 1981 (1981). 

[8] X.-L. Qiu and P. Tong, "Large-scale velocity structures in 
turbulent thermal convection", Phys. Rev. E 64, 036304 
(2001). 

[9] M. Sano, X.-Z. Wu, and A. Libchaber, "Turbulence 
in helium-gas free convection", Phys. Rev. A 40, 6421 
(1989). 

[10] C. Sun, K.-Q. Xia, and P. Tong, "Three-dimensional flow 
structures and dynamics of turbulent thermal convection 
in a cylindrical cell", Phys. Rev. E 72, 026302 (2005). 

[11] Y. Tsuji, T. Mizuno, T. Mashiko, and M. Sano, "Mean 
wind in convective turbulence of Mercury", Phys. Rev. 



Lett. 94, 034501 (2005). 

[12] E. Brown, and G. Ahlers, "Rotations and cessations of 
the large-scale circulation in turbulent Rayleigh-Benard 
convection", J. Fluid Mech. in press (2006). 

[13] E. Brown, A. Nikolaenko, and G. Ahlers, "Reorienta- 
tion of the large-scale circulation in turbulent Rayleigh- 
Benard convection", Phys. Rev. Lett. 95, 084503 (2005). 

[14] J.E. Hart, S. Kittelman, and DR. Ohlsen, "Mean flow 
precession and temperature probability density functions 
in turbulent rotating convection," Phys. of Fluids 14, 955 
(2002). 

[15] A. A. Draad, and F.T.M Nieuwstadt, "The Earth's ro- 
tation and laminar pipe flow", J. Fluid Mech.361, 297 
(1998). 

[16] G. Ahlers, E. Brown, and A. Nikolaenko, "The search 
for slow transients, and the effect of imperfect vertical 
alignment, in turbulent Rayleigh-Benard convection. J. 
Fluid Mech 555 347 (2006). 

[17] E. Brown, A. Nikolaenko, D. Funfschilling, and G. 
Ahlers, "Heat transport in turbulent Rayleigh-Benard 
convection: Effect of finite top- and bottom-plate con- 
ductivity." Phys. Fluids 17^075108 (2005). 

[18] The caption of Fig. 4 of gives R = 2.9 x 10 9 . We 
assume that this is a misprint since the values of the 
Taylor and Rossby numbers cited there, together with the 
Prandtl number, imply that R ~ 2.9 x 10 , and since 
the authors stated previously in their paper that there 
was no rotation of the LSC for R as large as 9.2 x 10 9 . 

[19] E. Brown, D. Funfschilling, and G. Ahlers, "Reynolds 
numbers of the large-scale flow in turbulent Rayleigh- 
Benard convection", unpublished. 



17 



[20] S. Grossmann, and D. Lohse, "Prandtl and Rayleigh 
number dependence of the Reynolds number in turbulent 
thermal convection" , Phys. Rev. E 66, 01635 (2002). 

[21] V. Ambegaokar, and B.I. Halperin, "Voltage due to ther- 
mal noise in the dc Josephson effect", Phys. Rev. Lett. 
22 25 (1969). 

[22] There appears to be a typo in Eq. 9 of Ref. HJ]. The 
limits of integration on the last integral over 8' should 
read from 6 to 2n instead of to 2n. 

[23] H.A. Kramers, "Brownian motion in a field of force and 
the diffusion model of chemical reactions", Physica (Am- 
sterdam) 7 284 (1940). 



[24] H.D. Xi, Q. Zhou, and K.-Q. Xia, "Azimuthal motion of 
the mean wind in turbulent thermal convection", Sub- 
mitted to Phys. Rev. E . 

[25] C. Sun, H.-D. Xi, and K.-Q. Xia, "Azimuthal Symmetry, 
Flow Dynamics, and Heat transport in Turbulent Ther- 
mal Convection in a Cylinder with an Aspect Raio of 
0.5", Phys. Rev. Lett. 95, 074502 (2005). 

[26] K.-Q. Xia, C. Sun, S.-Q. Zhou, "Particle image velocime- 
try measurement of the velocity field in turbulent thermal 
convection", Phys. Rev. £ 68, 066303 (2003). 



